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Abstract 

This paper considers a new model of individual displacement, based on fish 
motion, the so-called Persistent Turning Walker (PTW) model, which involves an 
Ornstein-Uhlenbeck process on the curvature of the particle trajectory. The goal 
is to show that its large time and space scale dynamics is of diffusive type, and 
to provide an analytic expression of the diffusion coefficient. Two methods are 
investigated. In the first one, we compute the large time asymptotics of the variance 
of the individual stochastic trajectories. The second method is based on a diffusion 
approximation of the kinetic formulation of these stochastic trajectories. The kinetic 
model is a Fokker-Planck type equation posed in an extended phase-space involving 
the curvature among the kinetic variables. We show that both methods lead to the 
same value of the diffusion constant. We present some numerical simulations to 
illustrate the theoretical results. 
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1 Introduction 



This paper considers a new model of individual displacement, the so-called 'Persistent 
Turning Walker' (PTW) model, which has recently been introduced to describe fish be- 
havior [27J. The fish evolves with a velocity of constant magnitude and its trajectory 
is subject to random turns (i.e. random changes of curvature) on the one hand and to 
curvature relaxation to zero on the other hand. The random changes of curvature can 
be interpreted as a way for the fish to explore its surroundings while relaxation to zero 
curvature just expresses that the fish cannot sustain too strongly curved trajectories and 
when the curvature becomes too large, it gets back to a straight line trajectory. The com- 
bination of these two antagonist behaviors gives rise to an Ornstein-Uhlenbeck process on 
the curvature. The curvature is the time derivative of the director of the velocity, while 
the velocity itself is the time derivative of position. The PTW process collects all these 
considerations into a system of stochastic differential equations. 

This model is, to the knowledge of the authors, original, and has appeared for the 
first time in the works by Gautrais, Theraulaz, and coworkers [27j- The present paper 
considers the large time and space scale dynamics of a two-dimensional particle subject to 
this PTW process. It rigorously shows (in the mathematical sense) that, at large scales, 
the dynamics of the particle can be described by a diffusion process and it provides a 
formula for the diffusion coefficient. To prove this result, two methods are considered. 

In the first method, the stochastic differential system itself is considered and the 
variance of the position is shown to behave, at large times, like a linear function of time. 
The diffusion coefficient is identified as the slope of this linear function. Because the 
curvature and the velocity angle can be explicitly computed, an explicit formula for the 
diffusion coefficient, involving some special functions, can be obtained. 

The second method considers the forward Kolmogorov equation of the stochastic pro- 
cess. This equation gives the evolution of the probability distribution function of the 
particle in the extended phase space (position, velocity angle, curvature) as a function 
of time. It is a Fokker-Planck type equation. The passage from the microscopic to the 
macroscopic scales relies on a rescaling of the Kolmogorov equation. This rescaling de- 
pends on a small parameter e <C 1, which describes the ratio of the typical microscopic 
to macroscopic space units. After this rescaling, the problem has the typical form of the 
diffusion approximation of a kinetic problem (see references below). The goal is then to 
study the behaviour of the solution as e — 0. It is shown that the solution converges 
to some 'thermodynamical equilibrium' which is a Gaussian distribution of the curvature 
and a uniform distribution of the velocity angle. The equilibrium depends parametrically 
on the density which satisfies a spatial diffusion equation. 

Finally, the connection between the two methods is made by showing that the diffu- 
sion tensor in the second approach can be represented by a formula involving the solution 
of the stochastic differential equation of the first approach. Additionally, this represen- 
tation leads to explicit computations which show that the two formulas for the diffusion 
coefficient actually coincide. This seemingly innocuous result is actually quite powerful. 
Indeed, the diffusion approximation method leads to a non-explicit expression of the dif- 
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fusion coefficient, involving the moments of a particular solution of a stationary equation 
involving the leading order operator of the Fokker-Planck equation. That this non-explicit 
formula is equivalent to the explicit formula given by the stochastic trajectory method 
is by far not obvious. In this respect, the stochastic trajectory method is more powerful 
than the diffusion approximation approach, because it directly leads to the most simple 
expression of the diffusion constant. 

A third route could have been taken and has been dismissed. This third method would 
actually use the stochastic differential equation itself to perform the diffusion approxima- 
tion in the forward Kolmogorov equation. We have preferred to use partial differential 
equation techniques. One reason for this choice is that these techniques can be more 
easily extended to more complex situations. One typical example of these more complex 
situations are the nonlinear systems which are obtained when interactions between indi- 
vidual are included. The inclusion of interactions between individuals within the PTW 
model is actually work in progress. 

From the biological viewpoint, one should not restrict the content of the paper to the 
sole expression of the diffusion coefficient. Indeed, once interactions between individuals 
will be included in the PTW model, it is not clear at all that the explicit computations 
which led to this expression will still be tractable. In the absence of an explicit solution 
of the stochastic differential system, there is little grasp to get information about the 
large scale behaviour of the system. By contrast, the diffusion approximation approach 
gives a systematic tool to study the large scale behavior of such systems, in all kinds of 
situations, be they linear or nonlinear. By its flexibility and its versatility, the diffusion 
approximation approach is the method of choice to study these problems. 

One of the most popular models to describe fish behavior is the discrete Couzin- 
Vicsek algorithm (CVA) [H [El [301 119] (see also [21 [IHl [251 [371 SQl ST] for related models). 
For a large scale modeling of fish behavior, it is efficient to look at continuum models, 
which use macroscopic variables such as mean density, mean velocity and so on. Several 
such models based on phenomenological observations, exist (see e.g. [2S1 I3S1 HZl US])- 
Several attempts to derive continuum models from the CVA model are also reported in 
the literature [Ml [HI [15]. In [211 122], a derivation of a continuum model from a kinetic 
version of the CVA model is proposed. However, few Individual Based Models for fish have 
been validated against experimental data with a comparable care as in [27] for the PTW 
process. As such, the continuum model derived in this paper has a &m experimental 
basis, although further work needs certainly to be done to fully validate its biological 
foundations. Additional references on swarm aggregation and fish schooling can be found 
in [llj. Among other types of animal societies, ants have been the subject of numerous 
studies and the reader can refer (among others) to [331 [IS] , and references therein. 

The derivation of macroscopic models from particle behavior has been initiated by the 
seminal works of Boltzmann, and later Hilbert, Chapman and Enskog. We refer to [13] 
for a mathematical perspective and to flBj for an introduction to the subject from a mod- 
eling perspective. More recently, the derivation of macroscopic models from microscopic 
behavior has been very productive in other context like traffic [H [32] or supply-chains 
[3]. Diffusion approximation problems for kinetic equations have been widely studied 
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in the literature, in the context of neutron transport (see e.g. P, E]), semiconductors 
[a Eal EH Ha], plasmas [HI [191 ED] or polymeric fluids [18]. 

The outline of the paper is as follows : in section [21 the PTW process is introduced 
and the main results are stated. In section [3] the diffusion coefficient is obtained by di- 
rect inspection of the trajectories of the stochastic differential system. In section [D, the 
diffusion approximation of the forward Kolmogorov equation of the stochastic process is 
performed. Section [5] is devoted to proving that the trajectory method and the diffusion 
approximation method give rise to the same value of the diffusion coefficient. In section 
[HI the theoretical results are illustrated by and complemented with some numerical simu- 
lations. A conclusion is drawn in section [3 Several proofs of auxiliary results, which are 
inessential for the main discussion are collected in three appendices (A, B and C). 

2 The Persistent Turning Walker model: presenta- 
tion and main results 

The starting point of the present work is a new model of fish motion based on experimental 
data taken from experiments run in La Reunion islands during years 2001 and 2002 |27j . 
The studied species is a pelagic fish named Kuhlia Mugil. Its typical size ranges between 
20 and 25 cm. The first experiments have been made with a single fish in a basin of 4 
meters diameter during two minutes. A video records the positions of the fish every 12-th 
of a second (see figure [T]). Then the data have been statistically analyzed and a model 
has been extracted |27] . 
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Figure 1: One trajectory of a Kuhlia Mugil fish 
The conclusion of the statistical analysis is that the trajectories are well described by 
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the following systems of stochastic differential equations: 



dx 

de 

H 
dn 



cm, 



—an dt + b dB^ 



(2.1) 



(2.2) 
(2.3) 



where x = {xi,X2) G is the (two-dimensional) position vector of the (centroid of the) 
fish, T{d) = (cos 6 , sin 6) is the director of the velocity vector with the angle 6 & U = 
M/27rZ measured from the Xi direction, k G M is the curvature of the trajectory and dBt 
is the standard Brownian motion. The magnitude of the velocity is constant and denoted 
by c > 0. The constant a is a relaxation frequency and b quantifies the intensity of the 
random curvature jumps, b has the dimension of 1/{L\/T) where L and T stand for the 
dimensions of length and time. 

The K-dynamics is a standard Ornstein-Uhlenbeck process. The term "6 dBt" models a 
diffusion process in curvature space while the term "—a K dt" expresses the tendency of the 
individual to return to a straight line trajectory. The curvature cannot increase endlessly 
as a consequence of the diffusion process, but rather, must relax to zero and the relaxation 
is stronger as the curvature gets larger. This model has been called the Persistent Turning 
Walker model (PTW) because it allows large excursions of the curvature towards positive 
or negative values, during which the spinning of the trajectory persists for a certain time. 

We stress the difference with more standard diffusion processes (such as those suffered 
by photons in a diffusive medium), in which the Brownian motion acts on the velocity 
itself (or, in the case of a velocity of constant magnitude, on the angle 6). In this case, the 
diffusion process acts on the second derivative of the particle positions, and the associated 
kinetic equation is of Fokker-Planck type. This model of photon diffusion is also relevant 
for a certain number of animal species |39j . 

In the PTW model, the diffusion process acts on the curvature, i.e. on the third 
derivative of the position vector. An intuitive justification of the relevance of this model 
for animal behaviour is by considering the non-differentiability of the Brownian motion. 
Because of this feature, the photon diffusion process involves infinite second derivatives of 
the position, i.e. infinite forces. However, an animal body can only exert finite forces and 
the muscles act only in such a way that the velocity angle undergoes smooth variations. 
The PTW model precisely presents this feature of having smooth second order derivatives, 
i.e. smooth forces. 

Our goal in the present work is to study the large-scale dynamics of the stochastic 
differential system fl2.ip - fl2.3l) . This is best done in scaled variables, where the dimension- 
less parameters of the model are highlighted. We use to = as time unit, xq = ca~^ 
as space unit, and curvature unit, and we introduce the dimensionless time, 

space and curvature as t' = t/t^, x' = x/xq and k' = k/kq. For simplicity, we omit the 
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primes. In scaled variables, the PTW model is written: 



I ^ (2.5) 

rfK = -Kdt + V2adBt, (2.6) 

where the only dimensionless parameter left is a such that 

The meaning of is the following: hj \fa is the amplitude of a curvature change during 
a relaxation time a~^, while c/a is obviously the distance travelled by the particle during 
this time. The product of these two quantities is dimensionless and is equal to ^J2ci. It 
quantifies the strength of the curvature jumps relative to the other phenomena. 

The individual dynamics can be translated in terms of a probability distribution 
f{t,x,9,K,)dxd9dK, of finding particles at times t with position in small neighborhoods 
dxdOdn of position x, velocity angle 6 and curvature k. The link between the individ- 
ual dynamics and the evolution of the probability distribution / is given by the forward 
Kolmogorov equation : 

dtf + r- V,f + ^def - d^{nf) - a^d,2f = 0. (2.8) 

This equation is an exact transcription of the individual djTiamics, where the initial 
value /o at time t = is given by the probability distribution of the initial conditions of 
the stochastic differential system fl2.4l) - fl2.6p . For more detailed considerations about the 
forward Kolmogorov equation and its link with stochastic differential systems, we refer 
the reader to [38l [5] . 

In order to capture the macroscopic dynamics, two possible routes can be taken, using 
either the stochastic differential system fl2.4p - fl2.6l) or the partial differential equation 
( 12. 8p . In this work, we follow both routes and verify that they lead to the same large- 
scale behaviour. The advantage of working directly on the stochastic system is that it 
is simpler and it leads to explicit formulas. However, as soon as the system gets more 
complicated, and in particular nonlinear, explicit solutions can no longer be found and 
this methodology can hardly be pursued. On the other hand, the PDE approach, which, 
in the present case is more complicated, is also more systematic and more general. In 
particular, it is generally usable in the more complex nonlinear cases (see e.g. [211 122]). 
A particular important complex situation is the case of many interacting fish. In future 
work, we plan to extend the PTW model to populations of interacting fish and to use the 
PDE approach to extract the large-scale dynamics of the system. 

From the analysis of the individual trajectories, explicit exact expressions for k and 
6 in terms of stochastic integrals can be found. Unfortunately, there is no such explicit 
result for the position x{t), but we can calculate the first two moments of the probability 
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distribution of x{t) explicitly, using the expressions of k, and 6. We show that the mean of 
the position vector stays at the origin: K{x(t)} = (0, 0) (where E denotes the expectation 
over all sources of randomness, in the initial data and in the stochastic process) and that 
the variance grows asymptotically linearly in time. More exactly, we prove: 

Theorem 2.1 Under assumptions on the initial conditions that will be specified later on 
(see 113. 1\) - I!i3.4^ ), the solution of system ((^.^[ )-( fl7^) satisfies: 

POO 

Var{f(t)} 2Pt, with V= exp {-a'^{-l + s + e-')) ds. (2.9) 

Jo 

The notation Var is for the variance over all sources of randomness. The asymptotic 
linear growth of the variance (12. 9p suggests that the dynamics of the system is of diffusive 
type at large times with diffusion coefficient T). We can find an expression of T> in terms 
of special functions. Indeed, we have 

Proposition 2.2 The following expression holds true: 

9 

e 



V 



)" 7(a',a'), (2.10) 



where ^{z^u) is the incomplete gamma function: 

7(2, m) = / e~H'-^dt. (2.11) 



V has the following series representation: 

00 

n\ {n + a^) 



P = e"^V A (2.12) 

n=0 



To investigate the large scale dynamics of the solution of the kinetic equation (12.81) 
(the existence of which can be easily proved, see proposition 14.21) . we need to rescale the 
variables to the macroscopic scale. Indeed, in eq. (12. 8p . all the coefficients are supposed 
to be of order unity. This means that the time and space scales of the experiment are of 
the same order as the typical time and length scales involved in the dynamics, such as, 
the relaxation time or the inverse of the typical random curvature excursions. Of course, 
in most experiments, this is not true, since the duration of the experiment and the size 
of the experimental region are large compared with the time and length scales involved 
in the dynamics. 

To translate this observation, we change the space unit xq to a new space space unit 
Xq = Xo/e, where e -C 1 is a small parameter. This induces a change of variables x' = ex. 
We make a similar operation on the time unit t'^ = to/r], t' = r]t with 77 <^ 1. Now, the 
question of linking 77 to e is a subtle one and is largely determined by the nature of the 
asymptotic regime which is achieved by the system. In the present case, we expect that 
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the asymptotic regime will be of diffusive nature, in view of theorem 12.11 and so, we will 
investigate the so-called 'diffusion approximation' which involves a quadratic relationship 
between rj and e: t] = e^. 

For this reason, we introduce the diffusive rescaling: 

t' = eH ; f = ex, (2.13) 

and we make the following change of variable in the distribution / : 

The scaling of the magnitude of the distribution function is unnecessary, since the prob- 
lem is linear. However, it is chosen in order to preserve the total number of particles. 
Introducing (12.131) into (12. 8p leads to the following problem for f^: 

edtf + r • V^r + \[^^^er - d^i^Kf) - a^d^.f] = (2.14) 

In order to analyze the large-scale dynamics of (I2.14p . we need to investigate the limit 
£ — > 0. We show that converges to an equilibrium distribution function (i.e. a function 
which cancels the 0{e~^) term of (12.141) ) /° which depends parametrically on the particle 
density n°(x, t) and evolves according to a diffusion equation. More precisely, we prove: 



Theorem 2.3 Under hypothesis \4jJ\ on the initial data to be precised below, the solution 
of Ili2.14\ l converge weakly in a Banach space also to be specified below, (see ^.21\) ) X: 

^-^^ n° ^^'^^ in X weak star, (2.15) 
27r 

where M is a Gaussian distribution of the curvature with zero mean and variance o? (see 
4-41 ) o-nd nP = n^{x,t) is the solution of the system: 

9*71° + ■ J° = 0, (2.16) 
jO = -DVsn^, (2.17) 

where the initial datum Uq and the diffusion tensor D will be defined later on (see ( [^.c/il ) 
and lli4-28^ respectively) . 

The following theorem connects the two methods by showing that the tensor D is 
related to V given by (12.91) : 

Theorem 2.4 The tensor D defined by l[4-28 ) satisfies : 

D = jld, (2.18) 

where T> is given by Ii2.9\) and Id denotes the 2x2 identity tensor. 

This theorem confirms that the trajectory method and the asymptotic PDE method 
are equivalent. The factor 2 between the two coefficients comes from the dimension of 
the problem. Indeed, "D is the variance of \x\'^ = -|- P while D is the diffusion 
coefficient in a given direction. 
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3 Large-scale dynamics of the PTW model by the 
trajectory method 

In this section, we want to show theorem 12.11 and proposition 12. 2[ We first specify the 
initial conditions. First, we fix the starting point of the particle at the origin : 

f(t = 0) = (0,0). (3.1) 

We suppose that the initial velocity angle is uniformly distributed on the one-dimensional 
sphere, i.e. : 

rlB 

dP{ei=, = e} = -. (3.2) 

For the curvature, we make the following observation: eq. (12. 6p predicts that the process 
K{t) converges exponentially fast to its stationary state, which is a Gaussian distribution 
with zero mean and variance equal to ^38j- We denote such a Gaussian distribution by 
A/'(0,a^). For this reason, we suppose: 

dP{K\t=Q = k] =U{<d,a^){K). (3.3) 

The last hypothesis on the initial conditions is the following: 

The processes 6{t = 0), K{t = 0) and Bt are independents. (3.4) 

We stress that this choice of initial conditions is for simplicity only. Completely arbi- 
trary initial conditions would lead to the same large time behaviour, but the computation 
would be slightly more complicated. Since we are mainly interested in the explicit expres- 
sion of V, a choice of initial conditions which simplifies the calculations is legitimate. 

We begin by proving the following proposition: 

Proposition 3.1 The solution of the stochastic differential equation ( [^.^p -( f^7^) with ini- 
tial condition given by I{3.1\) - ^i3^ satisfies : 

E{f(t)} = (0,0), Vt>0, (3.5) 
Var{f(t)} = 2! it- s)ex^{-a^ [-1 + s + e'')) ds. (3.6) 

To prove this proposition, we first establish explicit formulae for the solutions of (12.51) 
and (12. 6p . The proof is deferred to appendix A. 



Lemma 3.2 The solution of the stochastic differential system 1(2. 5\] . Ii2.6\) with initial 
conditions liS. ^) - (3^\ ) is given by: 

9{t) = Bo + Ko - K{t) + V2aBt, (3.7) 

K{t) = e-^Ko + V2ae-^ / e" dB^. (3.8) 

Jo 
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Additionally, 

e{t) = e, + Kl (3.9) 

where Kq is a Gaussian random variable independent of 6q with zero mean and variance 
given by : 

(3^ = Var{ir*} = 2a\-l + t + e''). (3.10) 

Proof of proposition I3.lt Using Lemma 13.21 we can compute the first two moments 
of x(t). Let us start witli tlie computation of ttie mean. If we write x{t) = {xi{t) , X2(t)), 
we liave : ^ ^ 

^i(^) = / cos6{s)ds , X2{t) = / sm9{s)ds, 



JO JO 

and, computing tlie mean : 

E{xi{t)} = ^ 1^ cose{s) ds^ = E{cos^(s)} ds. 

Now, we can develop 6{s) using (13. 9p : 

E{cos^(s)} = E{cos(^o + ^o)} =E{cos^o 

cos Kq — sin 9q sin Kq } . 

By the independence of 9q and Kq we finally have: 

E{cose(s)} = E{cos^o}E{cosis:o'} - E{sin(^o)}E{sin iTg"} = 0, 

since the expectations of cos 6q and sin Oq over the uniform probability distribution on Oq 
are zero. Finally, we have E{xi(t)} = 0, and similarly for X2- This proves (13. 5p . 
Now for the variance of x(t), we write: 

Var{x(t)} = ¥.{xl{t) + xl{t)} = 2E{x^(t)}. (3.11) 

by the isotropy of the problem. Then, 



E{xl{t)} = e\ (^J^ cose {s)ds 



t rt 



^0 

ft 



E{cos 6'(s) cos 6(u)} dsdu 

2 du ds E{cose{s) cos e{u)}. 
Jo JO 

Since m > s, we can write 9{u) as follows : 

e{u) = ^0 + y ^^{z) dz + J k{z) dz = 9o + K^ + K^ 

10 



where Kq and K"^ are Gaussian random variables independent of 6q with zero mean and 
variances Pi and respectively, thanks to (13. 101) . Then, using standard identities for 
trigonometric functions, we get 

E{cos^(s) cos^(m)} = E{cos(0o + ^o) cos(^o + + ^s")} 

= \{cos{29, + 2K^, + K:)+cos{-K:)). 

But since Oq is independent of and K"^ we have E{cos(26'o + IK^ + K")} = since the 
mean of a cos(^o + C) over the uniform distribution of is zero whatever the value of C. 
Then : 

E{cos(^(s))cos(^(m))} = ^E{cos(-K,")} 



If 1 w 

- / cos{y)—= e '^'^^-^ dy 

2 



1 _lfl2 

-e 2 



Indeed, an elementary computation shows that for any Gaussian random variable Z with 
zero mean and variance a^, one has 

E{cos(Z)} = exp(-aV2).) (3.12) 

Thus, 

pt ru 

¥.{xi{tf} = exp (-a^ (-1 + |m- s| +e"l"""l)) c/sdM. 

Ju=0 Js=0 

Using the change of unknowns w = u — s and y = u and inverting the order of integration 
we find : ^ 

E{xi(t)2} = / it-w) exp {-a^ {-I + w + e""")) dw. 

Jw=0 

Using (13.111) . we finally find (13.61) . which ends the proof of the proposition. ■ 

In order to prove [2?n we investigate the behavior of the variance Var{a;(t)} (given by 
( 13.61) ') when t +oo. 



End of proof of Theorem I2.lt we write, thanks to (13.61) : 

Var{f(t)} - 2Vt =-2 se-"'(-^+'+""') ds - 2 f^-'^H-i+s+e^) ds. 

J s=0 J s=t 

We have to show that the difference is bounded independently of t. For the first term, we 
have: 



t 



se 

s=0 



't 

< I e^'se-^'^ds, 
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and integrating by parts, we find 



se 



s=0 



< 



-te-^' - V + A 



For the second term, we have: 



< t I e°'e-°'"(is < te"'- 



< Co. 



This proves that the difference is Var{a;(t)} — 2Vt is bounded independently of t and 
completes the proof. ■ 



We now prove Proposition 12.21 which gives an explicit approximation of the diffusion 
coefficient. This approximation is useful for practical simulations. 

Proof of Proposition [2T2t the change of variables t = a^e"** in the integral (12.91) leads to 
(12.101) . The series representation fl2.12p follows from a similar series representation of the 
incomplete gamma function (see e.g. formula (8.354) of [21] )■ The series representation 
can also be found by expanding the exponential in the integral (12. lip in power series. 
This point is left to the reader. ■ 



4 Large-scale dynamics of the PTW model through 
the diffusion approximation of the associated ki- 
netic equation 

4.1 Formal asymptotics 

In this section, for the reader's convenience, we give a formal proof of theorem 12.31 We 
write (I2.14P as follows: 

edtr + T-v^r + ^Ar = o (4.i) 

where we define the operator A acting on functions u{6, k) as follows: 

Au = kOqu — di^^Ku) — a'^di^2U, (4.2) 

The formal investigation of the limit e — usually starts by considering the Hilbert 
expansion (see e.g. [16] for the general theory or |20j for an application in the context of 
Fokker-Planck equations): 

r = f + ef' + Oie'), (4.3) 

with being independent of e and inserting it into (14.11) . Then, collecting all the the 
terms of comparable orders with respect to e, we are led to a sequence of equations. The 
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first one, corresponding to the leading 0{e^^) term is Af^ = 0, which means that /° lies 
in the kernel of A. In section 14.31 we show that the kernel of A is composed of functions 
of the form /°()f:, x, 6, k) = n°(t, x)M{k,) / (27r) where M{k) is a normalized Gaussian with 
zero mean and variance a^: 

M(«:) = ^=L-e-^^, (4.4) 



and nP{t, x) is a function still to be determined. 

In order to determine n°, we first integrate (14. ip with respect to [6, k) G 11 x M and 
use that J AudOdn = 0. Defining the density n'^{t,x) and the flux J'^{t,x) by 

n'{t,x)= [ f dnde, J'{t,x)= [ —T{e)dKd9. (4.5) 
J9,K Je,K ^ 

we find: 

dtn" + Vs -r = 0. (4.6) 

We note that this continuity equation is valid for all values of e. Then, letting e — > 0, we 
formally have ^ n^. If we prove that J° given by (12.171) is the limit of J^, as e ^ 0, 
then, we can pass to the limit in (14. 6 p and find (12.160 . 

System (I2.16P and (12.171) is a diffusion system, which completely determines n^{t,x), 
given its initial datum UqIx). Here, for simplicity, we assume that the initial datum for 
(14. ip is of the form f^{0, x, 6, k) = no(x)M(K)/(2'7r) and the resulting initial condition for 
n° is therefore rig = no (in this formal convergence proof, we admit the functions and the 
convergences are as smooth as required). 

So, the only points left in the proof are the existence of a limit for and the validity 
of (12.171) for J°. Note that the existence of a limit is not obvious because of the factor e 
at the denominator of the integral (14. 5 p defining J^. To prove that the limit exists, we 
use the Hilbert expansion (14. 3 p again and compute Since A is linear, collecting the 
terms of order 0(£°) leads to: 

-A/' = ^Mf.v,„». (4.7) 

/vr 

Again using the linearity of A and the fact that it operates only with respect to the {6, k) 
variables, we can write the solution of (14. 7p as = — x ■ V^jn^, where x = (Xi;X2) is a 
solution of the problem 

M(k) ^ , , 

Ax = ^r. (4.8) 

This equation must be understood componentwise (i.e xi is associated with ti = cos 6 
and X2 with T2 = sin 6'). Since the right-hand side of (14. 8 p has zero average with respect 
to {0,k), proposition 14.51 below shows that it has a unique solution, up to an element of 
the kernel of A. We can single out a unique solution by requesting that x has zero average 
with respect to {9, n) as well. Then, all solutions to (14. 7p can be written as 

f = -X- V.^n" + n\t, x)M{k) /(2ti) , (4.9) 
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where the second term of (14. 9p is an arbitrary element of the kernel of A. We shall see 
that the determination of is unnecessary. 

Now, inserting the Hilbert expansion fl4.3l) into the integral fl4.5p defining J^, we find: 

r{t,x) = - [ fT{e)dKde+ [ ff{e)dKde + o{e) 

^ Je,K Je,K 

+ / ff{e)dKd9 + 0{e), (4.10) 

Je,K 

because /° is independent of 9 and j t{6) dO = 0. Therefore, has a limit when e —>■ 
and this limit is given by 

J\t,x) = I f^T{9)dKd9, (4.11) 

To compute J" we insert expression 04.91) into 04.111) and find 

J0(t,x)=l {-x-Vgn° + n^M/{27i))T{9)dKd9, (4.12) 

The second term vanishes and the first one can be written 

j^{t,x) = - (^j^ f®xd9dt^ Vgn^, (4.13) 

which is nothing but formula 02.17P with the diffusivity tensor D given by 04.281) . 

This shows the formal convergence of the solution of the Fokker-Planck equation 02.140 
to that of the diffusion system fl216D . ^TTl) . 

Now, to make this proof rigorous, we need to justify all the formal convergences. In 
the framework of the Hilbert expansion, this requires to work out the regularity of the 
various terms of the expansion. This is doable and actually leads to stronger convergences 
than the one we are going to prove, but this is a bit technical (see e.g. pQ]). 

What we are going to do instead is proving a convergence result in a weaker topology 
without using the Hilbert expansion technique. The method is close to the so-called 
moment method, which consists in integrating the equation against suitable test functions. 
This convergence proof is developed in section but before that, we state an existence 
result for the original Fokker-Planck equation 04.10 . 

4.2 Functional setting and existence result 

We define the differential operator D acting on smooth functions /(/t) by : 

Df = dM) + a''d^-f- (4.14) 

We state some properties of the proofs of which are easy and left to the reader. We 
recall that M{k) denotes the normalized Gaussian with zero mean and variance 04. 4p . 
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Proposition 4.1 Let f and g be smooth functions decreasing at infinity. The following 
identities hold true: 



Df = a' 
Dfg 
Dff 



d_ 

dn 
dn 
M 

dn 
M 

Df = 0^3ce 



M 



dnM 



—a 



—a 



M 



<9k 



M 

l_ 

M 



M 



dn 



fDg 



dK 
M' 



dK < 0, 



f = cM. 



(4.15) 
(4.16) 

(4.17) 
(4.18) 



The first identity translates the fact that M is the stationary measure of the Ornstein- 
Uhlenbeck process. The second one that D is formally self-adjoint with respect to the 
measure dn/M. The third one shows that D is dissipative. The same inequality holds with 
any non-decreasing function 77(/), indeed, / D{f) T]{f) dn < 0. If ?7 is the logarithm 
function, the corresponding quantity would be the relative entropy dissipation of / with 
respect to M. Entropy plays an important role in kinetic theory (see [13] for a review). 
Finally, the last quantity states that the kernel of D is one-dimensional and spanned by 
M. 

Proposition 14.11 shows that the natural norm associated with this operator has a 
weight and that the natural semi- norm is given by the right-hand side of fl4.17p . 
This motivates the introduction of the following functional spaces, endowed with their 
naturally associated Hilbert structures and norms: 

/dOdK 
\u{e,K)\^——<+oo}, (4.19) 

V = \ueH/ [ M dJ^)\Kd9<+oo\, (4.20) 

Ll = L^iRlH), X = L\[0,T]xRlV). (4.21) 

Identifying H with its dual, with have a Hilbertian triple V G H G V , where V is the 
dual of V and all injections are continuous. They are not compact because V does not 
bring any regularity with respect to 9. 

The existence proof follows closely the existence proof of [15] (see appendix A of this 
reference) and for this reason, is omitted (see also [20j). The proof relies on an existence 
theorem due to J. L. Lions [35j . 

Proposition 4.2 Let e > 0. We assume that /o belongs to L\j defined by lj^4-21\ ). Then 
there exists a unique solution to {2.14) with initial datum f^ in the class of functions 



Y defined by : 

Y={feX/dtf + e-'f ■ V,f + e'^Kdef E X'} . 
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Moreover, we have the inequality for any T > 0: 

2 

dKd9dxdt=\\r{0)\\l2^. (4.22) 

Estimate (14.221) is obtained via a Green and a trace formula for functions belonging 
to Y which can be deduced from the one proved in [T5] . 

4.3 Rigorous asymptotics 

We first study operator A given by (14.20 . i.e. Af = udef — Df and state some properties 
which will be proved in appendix B. We view A as an unbounded operator on the Hilbert 
space H with domain D{A) given by: 

D{A) = {u{e, K)eV /AueH}. 

Lemma 4.3 Operator A is maximal monotone. Moreover its kernel (or Null-space) is 
given by: 

Ker(A) = {cM, c G M}, (4.23) 

with M defined by ^.4^ . 



^ ^0 Jx£,K 



Lemma 4.4 The adjoint A* of A in H is given by A* f = —ndef — Df . It is a Maximal 
monotone operator with domain D{A*) = D{A) and Ker{A*) = Ker{A). 



Proposition 4.5 Let g E H . Then, there exists u G D{A) such that 

Au = g, (4.24) 
if and only if g satisfies the following solvability condition: 

I g{e,K)dedK = 0. (4.25) 

Moreover, the solution u is unique up to a constant times M . A unique solution can be 
singled out by presribing the condition 



u{9, k) dOdn = 0. (4.26) 



The same lemma applies to the equation A*u = g. 
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As an application of this lemma, let x be the solution of : 

Ax = m^. (4.27) 

with t(&) = {cosd , sin^). Since r has zero average over 6 and k, x is well-defined and 
unique thanks to Proposition I4.5[ Then, we define the tensor D by: 



D= t{9) (g) xd9dK. (4.28) 
Je,K 

Note that, since ^T{9)M{hi) dOdn = 0, it would not change the value of D to add any 
element of Ker(A) to x- 

Lemma 4.6 Let R denote the reflection operator u{6, n) Ru{6, n) = u{6, —k). Then, 
X* = Rx the unique solution (satisfying ^.2(j^ ) of 



A*X* = m ^, (4.29) 

and we have 

D = I f{9)®x*d9dK. (4.30) 

Proof: Obviously, D commutes with R: DR = RD while nde anticommutes with R: 
Kde{Ru) = -R{Kdgu). Therefore, RA = A*R. Since the right-hand side of (14.271) is 
invariant by R, applying R to both sides of fl4.27p leads to (14.291) . Then, the change of 
variables n' = —n in the integral at the right-hand side of (I4.30p shows that it is equal to 
D. I 

To study the limit e — > 0, we make the following hypothesis on the initial conditions. 

Hypothesis 4.1 We suppose that the initial condition /q is uniformly hounded in L\j 
and converges weakly in L\j to f^ as e —* 0. 

We can now prove theorem 12. 3[ The initial datum for the diffusion system (12. 16^ (12.171) 
will be shown to be: 

nO^Q = nl= f /°(f , 9, k) dKd9. (4.31) 
Proof of theorem 12. 3t By hypothesis 14.11 inequality (14.221) implies : 



2 



dKd9dxdt<C, (4.32) 
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with C independent of e. So (/^)e is a bounded sequence in L°°(0,T, L|^) and satisfies 



M 

X,tl,K 



dndedxdt <Ce\ (4.33) 



for any time interval T (by the diagonal process, we will eventually be able to take an 
increasing sequence of times T tending to infinity, so that the result will be valid on the 
whole interval t G (0,oo)). Therefore, there exists /° G L°°(0,T, L^^) and a subsequence, 
still denoted by satisfying : 

f f in L°°(0, T, LIj) weak star . 

Furthermore, with (14.331) . we deduce that /° = C{x,6,t)M{K). Then, letting e ^ in 
(12.141) . we get that Af^ = in the distributional sense. This implies that C{x,6,t) is 
independent of 6 and we can write 

f(t,x,e,K) = n%t,x)^^, (4.34) 

the quantity n°(t, x) = J /°(t, x, 9, k) dO dn being the density associated with /°. 

Our next task is to show that satisfies the diffusion model (12.161) . (I2.17P with initial 
condition (I4.3ip . We first note that is a week solution of (I2.14p with initial condition 
/q in the following sense: satisfies: 

r [ j;ef r. ^ ^ 1 ,,dKd9dx . f „ dndOdx , , 
/ / f{-edt^-T-Vi^ + -A*{^)) ^ dt = e I /oVt=o ^ • (4.35) 



for all test functions ip in the space C^([0,T) x ]R| x x M^) of twice continuously 
differentiable functions with compact support in [0, T) x R| x He x M.^. Again, the trace 
at t = has a meaning, thanks to a trace formula for functions in Y which is proven in 

We recall the definition of the fiux (14.50 . We prove that has a weak limit as 
e — >■ 0. To this aim, in the weak formulation (14.351) . we take as a test function = 
4>(t,x) ■ x*(^? ^) with x* the auxiliary function defined as the solution to (14.291) and is 
a smooth compactly supported vector test function of {x, t) . Although ip does not have 
a compact support, a standard truncation argument (which is omitted here) can be used 
to bypass this restriction. This allows us to write: 

f rr£/ o ^ \\/7 ^ PF^^ dndOdx . f ..7 ^dtidOdx 

Jo Jx,e,K e 2-n M J^n M 



Taking the limit e — > 0, we find : 



T 



lim / I r-<j)dxdt = 2n I / f{f-WM-Y*) dt 
£^0 




dudOdx 



dxdt. 
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where the exponent T denotes the transpose of a matrix. Using (14.301) and taking the hmit 
£ — >• shows that J"^ converges weakly (in the distributional sense) towards J° satisfying 



^ j_^J^-$dxdt = jn^V^-[{DYpjdxdt. 

This last equation is the weak form of eq. (12.171) . 

Finally, to prove (12.161) . we apply the weak formulation (14.351) to a test function of the 
form Lp = (j){t,x)M{K), where again, (f){x,t) is a scalar, smooth and compactly supported 
test function of {x,t) in x [0,T). This gives : 

~ / / {{edt + {t ■'Vg))(p) dndOdxdt = e / fQ(j)t=odKd6dx. 
Jo Jx,e,K Jx,e,K 

Dividing by e and taking the limit £ ^ 0, we get : 

- [ I {n^dtf + J° ■ V^0) dxdt = ! nl(j)t=Q dx, 

Jo Jx Jx 

where njj is defined by (I4.3ip . This last equation is exactly the weak formulation of 
equation (12.161) . with initial datum tiq. This concludes the proof. ■ 



5 Equivalence of the two methods 

In this section, we show that both methods lead to the same value of the diffusion coeffi- 
cient (theorem 12.41) . 

The first step is to show that we can approximate the solution of equation (14.241) by 
the solution of the associated evolution equation. More precisely, in appendix C, we prove 
the following lemma : 



Lemma 5.1 Let g in H satisfying lji4-^5\ ) and Uoo in -D(v4) be the solution of ^4-24 ) 



satisfying ( 4-26 ). Letuo G D{A) satisfying ( 4-26 ). Then, the solution u{t) of the evolution 
problem: 

dtu = -Au + g , Ut=o = uo, (5.1) 
weakly converges to u^o in H as t tends to oo. 

With this lemma we can explicitly calculate the tensor D and prove the theorem 12.41 : 
Proof of theorem 12. 4t Let x(^) be the solution of 

da = -Ax + m^. xit = o) = o. (5.2) 
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Thanks to Lemma [5. 11 weakly converges to x iii -f^ when t oo. It follows tht : 

x{t)®TdKde^^^ I x®fdKde. (5.3) 



Let us consider the first component of x(t), which we denote by u{t) and the inte- 
grals J^gu(t) cos 9 dudO and gu(t) sin 9 dudO. Because u satisfies f l5.2p . it admits the 
following representation (see the proof of Lemma [5.11) : 



u{t)= i'tJcosO^] ds, 



where Tt is the semi-group generated by —A (see [42j). With this expression, we evaluate 
the integral of u{t) against cos^ : 

/ u{t) cose dude = [ [ Ts (cose ^^^^ ds cos 6 dKdO 
Jk,9 JKfiJo V 27r y 



Kfi Jo 

* M(k) 

cos e — — T* (cos e) dn de ds. 



Jnfi 27r 

where T* is the adjoint operator of T in L^(^, k) generated by —A* , where 

Note that we are referring here to the adjoint in the standard sense and not in the 
weighted space H. This is why A* does not coincide with A* defined in Lemma I4.4[ The 
semi-group admits a probabilistic representation: for all regular functions /(6', n) 

T:{f){e, k) = E{f{e{t), KitWo = e,Ko = 

where {K{t),e(t)) is the solution of the stochastic differential equation (12.51) . (12.61) . Using 
this representation, we have : 



M {k)T* {cos e)dK = M{K)E{coses\eo = e,Ko = K}dK 

J K 

= Ejcos^sl^^o = e, kq = Z}, 



where Z is a random variable independent of Bf with density M. Using lemma 13.21 we 
have : 

E{cose,|^o = e,Ko = Z} = E{cos(^ + Ys)}, 

with Ys a Gaussian random variable with zero mean and variance (3^ given by ( I3.10p . 
Then: 

E{cos(e + Ys)} = E{cose cosF, - sin^ sinFj = cos^ EjcosFj, 

because the density of Yg is even and implies that E{sin Yg} = 0. Finally using (I3.12p . we 
have: 

M{k) T*{cose) dn = cose e"^ 



/ 



20 



Then, the first integral is given by: 

/■* f cos 6 _/3s 1 _i^s 

u{t) cos 9 dn d9 = / / cos9e~^d9ds= / -e ds. 

Jo Je 27r Jo 2 

We can proceed similarly to evaluate the integral of u(t) against sin()f:). This gives: 

r r cos 9 sin 9 

/ u{t)sm9dKd9= / / K{cosYs} d9ds = 0. 

Jk,0 Jo Je 27r 

It remains to evaluate the integrals involving the second component of vector which 
we denote by v{t). By the same method as for u{t), we get : 

f f /■* 1 _£i 

/ v(t) cos 9 dKd9 = and / v(t) sm9 dKd9 = / -e ^ dt. 

J Kfi J K,e Jo 2 

Collecting these formulae, we can write: 

x(t)®rdKd9 = '^^Id, 

with T>{t) = Jq Q-a^i-i+u+e ") Taking the limit t +oo and using equation fl5.3p . 
shows that (I2.18p holds true and completes the proof of the theorem. ■ 



6 Numerical simulation 

We simulate individual trajectories satisfying equation fl2.4p - fl2.6l) with initial conditions 
given by f l3.ip -( l3l4l) . If we fix a time step At, using (13 .70 . fl3.8p . we have : 

l^{n+l)At = 7^nAt + G{n+1) 

9(n+l)M = 9o + Kq — /«(n+l)Ai + V^"-B{n+l)At • 

with 7 = e~^* and G(^n+i) a Gaussian random variable with zero mean and variance 
2a^ (l — e~^^*) independent of K„At- With this formula, we can simulate recursively the 
process (fi;„At, 9nAt)n exactly (in the sense that it has the same law as the exact solution). 
To generate the Brownian motion, we just compute the increments -B(„+i)At — BnAt since 
they are Gaussian and independent of BnAt- On the other hand, these increments are 
not independent of G(^n+i)- Fortunately, we can compute the covariance matrix of the 
Gaussian vector (G'(„+i), 5(„+i)At - B^At) ■ 

ij(n+l)At — iJnAt J 
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where A/'(0, C) is a two-dimensional Gaussian vector with zero mean and covariance matrix 
C given by: 

[ V2a{l - e-^*) At 

Knowing this covariance matrix, we can simulate the Gaussian vector {G(^n+i), -B(n+i)At — 
BnAt) using the Cholesky method : we generate (Xi,X2) a vector of two independent 
normal law, and take ^/C{Xi, X2)^ as realization of the Gaussian vector. 

Now for the position x, since we do not have any explicit expression, we use a discrete 
approximation scheme of order 0((At)^). For example, the first component xi of x is 
approximated by: 

Xi{{n + l)At) = xi{nAt) + cose{s)ds 

JnAt 

fti xi{nAt) + ^{cose{nAt) + cos6'((n + l)At)). 

We present four trajectories obtained with different values of the parameter a in 
figure [21 As the parameter a increases, the excursions towards large positive or negative 
curvatures become larger. As a consequence, the spinning of the trajectory around itself 
increases and, from almost a straight line when a = 0.1, the trajectory shrinks and looks 
closer and closer to a wool ball. We can figure out that the diffusion coefficient decreases 
as a increases which is confirmed by formula (12. 9p since the function — l + s + exp(— s) > 
for s > 0. Additionally, it is easily seen that V ^ oo like when a —>■ while V 
when a ^ CO. 

To illustrate theorem 12.11 we use a Monte-Carlo method to simulate the variance of 
the process x. We simulate N independent trajectories and we compute the variance 
of the sample at each time step. In figure ([3]), we compare the result obtained with 
= 2000 and the theoretical prediction given by the (13. 6p . The figure shows an excellent 
agreement between the computation and the theoretical prediction. Additionally, after an 
initial transient, the growth of the variance is linear, in accordance with the theoretical 
result (12. 9p . We can use the slope of the asymptotically linear part of the curve to give a 
numerical estimate of the diffusion coefficient V. Fro this purpose, we fit a straight line 
(in the mean-square sense) between times T/2 and T. We remove the data between and 
T/2 because the initial transient is not linear and including them would deteriorate the 
accuracy of the measurement. We compare the slope of the fitted line with the theoretical 
value (12. 9p . We report the result of this comparison for two values of a (a = 0.1 and 
a = 2., with T = 1200 time units in Table [H The approximation is quite good, with an 
error comprised between 2 and 3%, which can be attributed to numerical noise and to an 
unsufficient approximatin of the asymptotic state. To illustrate the infiuence of the initial 
transient, we take T = 120 time units in the case a = 0.1 and report the result in Table 
121 There, the approximation is quite poor, because the asymptotic state has not yet been 
reached. 

In order to illustrate theorem l2.3[ we plot the spatial density n{t, x) of the distribution 
f{t, X, 6, k) using a Monte-Carlo algorithm for a = 2 and T = 30 time units on figure |H 
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a=0.1 a=0.5 




-20 -10 10 20 -20 -10 10 20 



Figure 2: Four trajectories simulated with different value of a: a = 0.1 (top left), a = 0.5 
(top right), a = 1. (bottom left) and a = 2. (bottom right). The simulation is run during 
120 time units with a time step dt = 0.05 time unit 





D simulation 


D theoretical 


relative error 


a = 0.1 


98.5 


101 


2.5 % 


a = 2 


0.708 


0.725 


2.4 % 



Table 1: Diffusion coefficient: comparison of the numerical estimate obtained by fitting 
the numerical values with a straight line over the time interval [T/2, T] with the theoretical 
prediction (12.91) . T = 1200 units of time. 





D simulation 


D theoretical 


relative error 


a = 0.1 


58.8 


101 


72 % 



Table 2: Same as Tabled] but with T = 120 units of time. The agreement is poor because 
the asymptotic state is not reached. 
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Variance of the process a=0.1 



Variance of the process a=2 
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Theoretical 




20 
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(a) 



100 120 



200 
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100 




40 60 80 100 
times in second 



(b) 



120 



Figure 3: Variance of the process x{t): comparison between the numerical simulation 
(points) and the theoretical prediction (solid line) for two values of a : a = 0.1 (left) and 
a = 2 (right). 



We see that the density has the Gaussian shape of the solution of a diffusion equation, in 
accordance with the prediction of the theorem. To make a more quantitative comparison, 




-30 -30 



Figure 4: The spatial distribution n for a = 2 at time T = 30 time units. 

we compare it with the asymptotic prediction, i.e. the solution of the diffusion equation 
(12.161) . fl2.17p . by computing the difference in norm. The results are reported in figure 
O We plot the norm of the difference for a = 1. and for four values of the parameter 
e : 1, 1 1 and ^. As expected, the agreement is better as e is smaller. However, at large 
times, all solutions are eventually close to the solution of the diffusion equation. Roughly 
speaking, the time at which the solution of the diffusion equation starts to be a good 
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approximation of the solution of the kinetic equation scales like e. This means that, after 
an initial transient the duration of which may depend on e, the solution is close to that 
of the diffusion equation, no matter the value of e. 



1.8 
1.6 

1.4 

1— 1.2 

_J 

O 

1 

I 

_E_ 0.8 
0.6 
0.4 
0.2 


5 10 15 20 25 30 35 

times in seconds 

Figure 5: norm of the difference between n^(t, ) = J /"^(t, ■, 6', k) dO dn and its asymp- 
totic limit 'nP{t) as a function of time t, with a = 1 and N = 10^ simulation particles. 
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7 Conclusion 



In this paper, the large-scale dynamics of the 'Persistent Turning Walker' (PTW) model 
of fich behavior has been analyzed. It has been shown, by two different methods, that the 
large scale limit of this model is of diffusion type, and an explicit formula for the diffusion 
coefficient has been provided. While the direct analysis of the stochastic trajectories 
provides a direct route to the value of the diffusion constant, the diffusion approximation 
of the associated forward Kolmogorov equation, which is of Fokker-Planck type, gives 
a more systematic way to extend the theory to more complex nonlinear cases. Such a 
nonlinear situation will be encountered when, in the near future, the nonlinear interactions 
between the individuals will be introduced within the PTW model. We expect that, in this 
context, the diffusion approximation methodology will have to be exploited thoroughly 
to allow access to the large scale behaviour of the system. 
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Appendix A: proofs of section [3 

Proof of Lemma 13. 2t formula (13. 8p is standard in the theory of Ornstein-Uhlenbeck 
processes [38] • To obtain (13.71) . we integrate ^(t) with respect to time: 

/* i ft/ I* s 

I k{s) ds = {l- e-*)Ko + V2a I I Q-'e"" dB^ds. 
Jo Jo Jo 

Interchanging the order of integrations and integrating with respect to s, we deduce: 

[ K{s)ds={l-e-')Ko + V2a [ (1 - e^^*"")) 
Jo Jo 

Then we develop the integral : 

/ k{s) ds = {1 - e"*)/to + V2aBt - (/t(t) - e~*Ko)- 
Jo 



This formula can be rewritten: 



/ k{s) ds = Kf) — K{t) + \/2aBt, 
Jo 

which easily leads to (13.71) . 

We now calculate the mean and the variance of Kq = k{s) ds. Since k{s) is of zero 
mean, its integral Kq is also of zero mean: K{Kq} = 0. Now for the variance of Kq, we 
can write : 

Var{ir*} = E I ^(1 - e-')Ko - v^ae"* ^ e' dB^ + V2aBt^ | . 
Using that kq and Bg are independent, we can develop the square and get: 

Var{K*} = (1 - e''YE{Kl} + 2a'^E | (^-e"* ^ e' dBs + Bt^ 
Let us consider the second term. By Ito's formula, we have 

E I (^-e-' e' dBs + ^t) | = e"^*E | (^j^ e' dB^j | - 2e-*E ^ e' c/5,| + E {B 

= e-2* ^ e^' ds - 2e-*E |st(e*fit - ^ e'B^ ds^+t, 

where the Ito correction term is zero due to the fact that exp s is a deterministic process. 
We can simplify this expression again since K{BtBs} = min(t, s) and get: 

E <; ( -e"* f e' dBs + B^\ = ^—^ 2e-* ( eH - / e's ds\ + t 

1 - e~^* 
2 
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2 (1 -e^*) +t. 



Using also that 1E{ko} — ^ ^ variance of K\ is written: 

Var{K*} = (1 - e-*)2«2 + (^^^ "2(1- e"*) + . 
Developing and simplifying the expression, we find fl3.10p . which ends the proof. 
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Appendix B: proofs of section 



4 



Proof of Lemma 14. 3t Let u G D{A). Then, ndgu G V and Lemma Al of [15] shows 
that the Green formula for functions u & V such that udeu G V is legitimate. Therefore, 
taking the inner product of A{u) against u, we find: 

^ dedK > 0. (7.1) 

So, A is a monotone operator in H. To show that A is maximal monotone, we prove that 
for any g & H, there exists u G D{A) such that : 

u + Au = g. (7.2) 

Taking the inner product of (17. 2p against a test function ip in the space ViJlQ x M^) 
of infinitely differentiable and compactly supported functions on Hg x Mf, leads to the 
variational problem : 

lju (, - .ft.) -i + Ma. (^) a. (^)] rf^d. = £ 9^ (7,3) 

Again, the same theory as in the appendix A of [15j (based the result by J. L. Lions in 
[35]) applies to prove the existence of a solution to (17. 3p with u in V such that ndgu G V. 
From there, it immediately follows that u G D{A). 

It is immediate to see that any function of the form u{6, k) = CM{k) for any constant 
C belongs to the kernel of A. Conversely, suppose that u G Ker A. Then, by (17. ip . there 
exists a function C{0) G LF'iJ^) such that u{9,k) = C{9)M{k). But again, A(u) = 
implies that KdgC{9) M = 0. So C{6) is a constant, which proves (14.230 . ■ 



<A{u),u>H= / a^M d^{^ 



u 
M 



Proof of Proposition 14. 5t the 'only if part of the theorem is obvious since, us- 
ing Green's formula (again, obtained by adapting that of appendix B of [15], we have 
/ Au d9 dK = 0. 

To prove the 'if part, we borrow a method from (for instance) [12]. To find a solution 
to (I4.24p . we look at a perturbed equation : 

Xu + Au = g, (7.4) 

with A > 0. Since A is maximal monotone in H (Lemma 14. 3p . eq. (17. 4p admits a solution 
u\ for all positive A ([9|). To prove the existence of a solution to ( I4.24p . we want to extract 
a subsequence, still denoted by {u\) which converges weakly in H. For this purpose, it is 
enough to show that there exists a bounded subsequence. 

We proceed by contradiction, supposing that the (full) sequence Nx = \\ux\\h +oo. 
We define U\ = U\ satisfies ||f/A||H = 1 for all A and 

XUx + AUx = 4r- (7-5) 
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Since {Ux)\ is a bounded sequence in H, we can extract a subsequence (still denoted by 
Ux) such that Ux ^ U in H weak as A — > 0. Taking the limit A — > in (17.51) . gives 
A{U) = 0. If we take the inner product of fl7.5p with Ux and then pass to the limit 
A — > 0, we also find that U belongs to V. So Lemma [4.31 applies and gives U = cM with 
a constant c G M. Using fl4.25p . we also have : 

< XUx + AUx, M >H=< ^, M >H= 0. 

So f aUxdOdn = for all A. Taking the limit A — > leads to f „U dOdn = 
J^gCM{K,) dOdn = 2ttC = 0, which implies U = 0. This proves: 



Ux 







in H weak 



(7.6) 



To get a contradiction, we now prove that the convergence is strong. 

To this aim, we introduce a decomposition of the space H into two orthogonal sub- 
spaces. Let L be the closed subspace of H defined by : 

L = {c{e)M/ c{e)eL\Ue)}, 



with M defined by (14.41) . So H = L ® L . We also define the orthogonal projector P 
of H onto L such that Pf = (Xc fi^^ ^) d'^) Using this projection, we decompose the 
sequence {Ux)x as follows: 

Ux = cx{e)M + Vx, (7.7) 



with Vx G L-^i i.e. J^vxdn = 0. To demonstrate that Ux — in if strongly, we first 

demonstrate that vx — ^ in ff strongly. 

Taking the inner product of the equation satisfied by Ux (17.51) with Ux gives : 



M\Ux\\l + 







I M 


'e,K 





d9dK = ^ < g,Ux>H ■ 
^Va 



Since d^j^ = ^^"^ l|f^A||// = 1; we get by taking the limit A ^ 



■M 



I M 




le,K 





dOdn 0. 



(7i 



Now Gross inequality [31j gives, for any v 

2 



<9k 



M 



MdK+ [ fdK > I/I 



Then, since vx dhz = 0, we deduce: 



di-i 
M' 



(7.9) 



vx_ 
M 



MdK> 
30 



\vx\' 
M 



dn. 



Integrating this inequality with respect to 6 and using (17. 8p . we find: 

II^aII// 0, in H strong. 



(7.10) 



To prove the convergence of ca, we define the bounded operator T: H ^ L'^(Jlg) such 
that Tf = nf dn. Having T acting on (17.5P and taking the hmit A — > 0, leads to: 



TAUx^Q in (n^) strong 
If we develop the left-hand side, we find: 

TAUx = 



(7.11) 



' K^deUxdn— / [K,dK{KUx) — a'^KdK2Ux] dn 

K J K 



K^deUx dn + kUx dn. 



But using the decomposition Ux = cxM + vx (17. 7p . we have: 



/ KUxdK 




/ KVx dK 


^0 


J K 


LH0) 


J K 


L\e) 



So, (ITTTil leads to : 



/ K^deUxdK ^ in L?'(Jle) strong. 

J K 



(7.12) 

If we define hx{0) = J^K^UxdK, (17.121) is equivalent to saying that ||(96)/;.a||l2(6») 0. 
Using the Poincare-Wirtinger inequality [9], there exists a constant Cq such that: 



— hx\\L2(0) < Co\\d0hx\\L2(e), 



(7.13) 



with hx = /o^'^^a(^) dO. Then, we develop hx- We get: 




A^O 



K^Ux dude =< Ux, Mk^ >h — ^ in R, 



Jk 



since Ux converges weakly to zero (see (17. 6p ). So, (I7.13P leads to hx ^— ^ in L'^{ll0) 
strong. If we develop hx we find: 

hx{9) = / k'^{cx{0)M + vx)dK = a'^cx{9) + / K^vxdn. 

J K J K 

Now, K^vxdK converges to zero in L'^{6) strong because of (I7.10p and we finally have : 



A^O 



cx{0) ^ in L^e) strong. 
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Using the convergence of ca and vx, we can now prove the strong convergence of Ux 
to in H: 

WUxWl = \\cxM\\l + WvxWl = ||c,||i.(,) + WvxWl ^ 0, 

which contradicts the fact that Ux has unit norm in H. This shows that there exists a 
bounded subsequence in the sequence ux- In fact, since the same proof can be apphed to 
any subsequence, this shows that the whole sequence ux is bounded, but this is useless 
for our purpose. 

We conclude the proof of Proposition 14.51 as follows: there exists a subsequence ux and 
a function m in H such that ux ^ u in H weak. Taking the limit of (17.41) as A —> 0, we 
deduce that Au = g in the sense of distributions. However, since g E H, eq. Au = g also 
holds in H. Moreover if we take the inner product of ( 17. 4p with ux and pass to the limit 
A — >• 0, we find that u belongs to V. So u belongs to D{A), which ends the proof of the 
'if part of the statement. 

Finally, to prove uniqueness, we just remark that, two solutions of (14.241) differ from 
an element of the kernel of A and we apply (14.231) . This ends the proof. ■ 
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Appendix C: proofs of section [5 



Proof of Lemma I5.lt the proof borrows some ideas from [23], but is simpler, due to 
the hnear character of the problem. The difficulty is getting some compactness in time. 
Here, instead of considering time translates of the solution as in p4], we will consider 
time integrals over a fixed interval length At. 

Since operator A is maximal monotone on H (see Lemma US]), operator —A generates 
a semi-group of contractions Tf on H. Moreover the solution of (15. ip is given by: 

u{t) = Tj(mo) + / Ts{g) ds. 
Jo 



We define f{t) = u{t) — u^o which satisfies : 

dtf = -Af, /i=o 



(7.1) 



with fo = Uq — Moo and /o(ft:) = 0. To prove the weak convergence of u(t) to u^o, we 
have to prove that f{t) converges to zero weakly in H. 

To this aim, we make an orthogonal decomposition of f{t) as in the proof of Proposition 
1431 f{t) = c{t)M + v{t), with c{t) e L'^ijle), v{t) E H and J^v{t) dn = 0. Taking the 
inner product of (17. ip with /, we get : 



n 2 



dndO. 



Using the decomposition of f{t) and noticing that (9^ (-^) = this equality be- 



comes: 



■dtiWc 



L2 



+ \\v 



I) 



vjtl 
M 



1 2 



dndO. 



(7.2) 



If we apply the Gross inequality (17.91) . we get: 



2^* (lie 



L2 



+ \\v 



i)<-\\vmH- 



Since c(t) is bounded by ||/o||//, by integrating with respect to time, we have : 



1„ 
2I"" 



< - / \\v{s)\\]jds + C. 



Using the Gronwall lemma, we deduce that v{t) decays exponentially fast to zero strongly 
in H: 

v{t) Q m. H strong. 



It remains to prove the convergence of c(t) to zero. We integrate (17. ip with respect to 
n. This gives, using that M(k) dn = 1 and v(t) dn = : 



(7.3) 



dtc{t) = 80 Kv{t) dn. 
J ft 
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Now if we pre-multiply by k before integrating with respect to k, we obtain : 

dt / Kv{t) dn = a^decit) + dg K^v{t) dn — Kv{t) dn. (7.4) 



We fix a time interval At and integrate fl7.4p over this time interval. This leads to: 

pt+At p pt+At 

K{v{t + At) — v{t)) dn = Q^de / c{s) ds + de v{s)dsdK 

J t J K J t 

t+At 

K I vis)dsdK. 



' K J t 

Since v{t) converges to zero in H , we have, in the sense of distributions: 

/t+At 
c{s) ds 0. (7.5) 

Since c belongs to L°°((0, oo)*, L'^iUe)) (see (E2D, we have //+^* c{s) ds G L°°((0, oo)*, L'^{Ue)). 
So there exists a subsequence such that f^~^^^ c(s) ds is weakly convergent in L'^ijle). Ac- 
tually, (17.51) implies that there exists a constant function with respect to 9, depending on 
At and denoted by L{At) such that 

rt+At 

c{s) rfs*^-i°'L(At). 

To deduce the convergence of c(t), we have to control the derivative of c{t) in time. For 
this purpose, we rewrite : 

"t+At pAt / \ 

c{s) ds = I c(t) + / dtc{t + z) dzj ds 

pAt PS p 

= Atc{t)+ / '^e Kv{t + z) dndzds. 

Jo Jo Jk 

Using again the convergence of v{t) to zero, we find : 

Atc(t) *^°°L(At), 

or defining the constant C = we have c(t) C in L'^ijlg) weak. 

To complete the proof, it remains to prove that C is equal to zero. Now, since eq. 
(17. ip is mass preserving i.e.: 



dt / fit) dKde = - / Af{t) dude = 0, 
Jk^ J K,e 

we have ^ f{t) dudO = /(O) dudO = 0. Also : 

f{t) dndO = [ {c{t)M + v{t)) dndO [ C d9 = 27rC. 



So C = 0. This proves f{t) in if weak and completes the proof. 
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